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ABSTRACT 



One possible scenario for the formation of carbon-enhanced metal-poor stars is the accretion of carbon-rich material from a binary 
companion which may no longer visible. It is generally assumed that the accreted material remains on the surface of the star and does 
not mix with the interior until first dredge-up. However, thermohaline mixing should mix the accreted material with the original stellar 
material as it has a higher mean molecular weight. We investigate the effect that this has on the surface abundances by modelling a 
binary system of metallicity Z = 10~ 4 with a 2M G primary star and a 0.74 M Q secondary star in an initial orbit of 4000 days. The 
accretion of material from the wind of the primary leads to the formation of a carbon-rich secondary. We find that the accreted material 
mixes fairly rapidly throughout 90% of the star, with important consequences for the surface composition. Models with thermohaline 
mixing predict very different surface abundances after first dredge-up compared to canonical models of stellar evolution. 
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1. Introduction 

Carbon-enhanced, metal-poor (CEMP) stars a re defined as stars 
with [C/Fe]>+1.0 dBeers & Christliebl I20051 A possible for- 
mation scenario for these stars is the accretion of carbon- 
rich material from a companion. It is c ommonly assumed (e.g. 
ICavrel et ail 12004 iLucatello et~aT1l2005l) that this accreted ma- 
terial sits unmixed on the surface of the star up to the point at 
which the star undergoes first dredge-up on the red giant branch. 
At this point, it is thought that the accreted material becomes 
mixed with part of the interior of the star in the deep convective 
envelope that develops. However, this picture neglects thermo- 
haline mixing. 

Thermohaline mixing is the process that occurs when the 
mean molecular weight of the stellar gas increases towards the 
surface. Although this situation is dynamically stable as long as 
the density decreases outward, a secular instability sets in when 
a gas element is displaced downwards and compressed. Such an 
element must be hotter than the surrounding gas by virtue of 
its higher molecular weight and if it loses heat to its surround- 
ings its density increases and it will continue to sink. As a result 
mixing occurs on a thermal ti mescale, until the mole cular weight 
difference has disappeared (Kippenhah n et al. l ll980b . In interact- 
ing binary systems, mass transfer (by wind accretion or Roche 
lobe overflow) may lead to the transfer of material that has been 
processed by nuclear burning from the primary to the surface of 
the less evolved secondary. This alters the surface composition 
of the secondary. The extent to which the surface composition is 
changed depends on how much the transferred matter is mixed 
with the pristine stellar material. 

We expect thermohaline mixing to occur in a star that has 
accreted material from a carbon-rich companion as this mate- 
rial will have a greater mean molecular weight than that of the 



unevolved star. The aim of this work is to determine the effect 
that thermohaline mixing has on the surface abundances of such 
an object. In section 2, we describe our model. In section 3, we 
describe the results and in section 4 these are discussed in the 
context of observations of CEMP stars. 



2. The model 
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We use the rapid synthetic binary code of llzzard et al.l d2006l) to 
model the evolution of a binary system of metallicity Z = 10~ 4 
consisting of a 2 M primary with a 0.74 M Q companion in a cir- 
cular orbit with an initial period of 4000 days. The secondary 
mass is chosen as being representative of a companion with a 
mass close to the turn-off mass of the halo. The primary mass 
was chosen as an example of an asymptotic giant branch (AGB) 
star that becomes very C-rich. During the thermally pulsing 
asymptotic giant branch (TP- AGB) the primary undergoes third 
dredge-up and becomes carbon rich while the en velope is eroded 
by th e stellar wind. Based on the Bondi-Hoyle (Bon di & Hovld 
1 19441) pr escription for wind accretion (as used in the synthetic 
code, see llzzard et al.ll2006l and references therein), we find that 
the secondary accretes 0.09 M of material from the wind of the 
primary and we record the average composition of this material 
(see Table [TJ. The total mass of material accreted and the av- 
erage composition of this material are then used in creating a 
detailed model of the secondary. 

We model the secondary using a version of the st ellar evo- 
lution code originally developed by Eggleto^ d!971l) and up- 
dated by many authors (e.g. Po ls et al] 1 19951). Our vers i on in- 
clud es the nucl e osynth esis subroutines of Stancliffe et al. (2005) 
and [Stancliffe (200 5|). T he code uses the opacity routines 
of lEldridge & Toutl d2004l). which emp loy interpolation in the 
OPAL tables (Iglesias & Rogers 1996) and which account for 
the variation in opacity as the C and O content of the material 
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varies. We create a 0.74 M model and evolve it from the pre- 
main sequence to the age of the system predicted by the rapid 
code. The initial abundances are cho s en to b e solar-scaled based 
on the values of lAnders & Grevessd d 19891) . with [Fe/Hl3=-2.3. 
A mixing length parameter of a = 2.0 was used and no convec- 
tive overshooting was employed. 

Thermo haline mixing was inclu ded according to the pre- 
scription of Kippen hahn et al.l dl980h . They model it as a diffu- 
sive process with a diffusion coefficient, Z) m , that can be written 
as 



Ah = 



l6acT 3 H P 

(V ad - V)Cpp* 



(1) 



where Hp is the pressure scale height, a is the radiation pressure 
constant, c is the speed of light, T is the temperature, V a d is the 
adiabatic temperature gradient, V is the actual temperature gradi- 
ent, cp is the specific heat at constant pressure, p is the density, k 
is the opacity, r is the radius and yiz is the mean molecular weight 
of the material. 

Once we have a model of the desired age (specifically 
8.6 x 10 8 yrs), we create two model sequences. In both, we ac- 
crete 0.09 M of material with the composition given by the 
rapid code. The accretion rate is set to lO _6 M yr _1 , which is 
approximately the rate at which the secondary would accrete ma- 
terial from the AGB superwind. In one sequence we include the 
effects of thermohaline mixing; in the other we do not. 



3. Results 

A Hertzsprung-Russell (HR) diagram showing the evolution of 
the two models is presented in FigureQ] The model without ther- 
mohaline mixing (dashed line) is cooler than the model with 
thermohaline mixing (solid line). This is because the accreted 
material is not mixed into the interior of the star. Although the 
star has a surface convection zone, this is too shallow to cause 
significant mixing. Because the accreted material has a higher 
opacity on account of it being extremely carbon rich, the surface 
temperature of the model without thermohaline mixing is lower 
and the star's position is shifted in the HR diagram relative to 
that of the model with thermohaline mixing. 

Figure [2] shows the hydrogen abundance profile of the star 
just after accretion for both models. Before accretion the surface 
material has a mean molecular weight of /u = 0.593 whilst the 
accreted material has a mean molecular weight of ft — 0.657 
(mainly caused by its higher helium abundance). Thermohaline 
mixing is strong enough to affect the profile during the accretion 
phase so the model with thermohaline mixing (dashed line) has 
a slightly shallower hydrogen abundance profile near the surface 
after accretion than the unmixed model (solid line). Figure|2]also 
shows the H-abundance profile when thermohaline mixing has 
ceased (dotted line). Note that by the time this state is reached, 
significant H-burning has occurred in the core. 

We find that thermohaline mixing results in the star being 
mixed down to within 0.1 M of the centre. This represents al- 
most 90% of the star. Only the very central material escapes be- 
ing mixed as it is here that material has experienced sufficient 
nuclear burning to raise its mean molecular weight. The mixing 
is rapid and the star reaches its equilibrium configuration within 
about 10 9 years after accretion. This is about one-tenth of its 
main-sequence lifetime. 



1 This is the logarithm of the ratio of the iron abundance to the hydro- 
gen abundance, relative to the solar value, i.e. [Fe/H] = \og(Xp c /Xn) - 
\og(X F JX H ) . 
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Fig. 1. Hertzsprung-Russell diagram showing the evolution of 
the two models. The model sequence with thermohaline mixing 
is the solid line; the dashed line is the model sequence without 
thermohaline mixing. The excursion at the beginning of the evo- 
lutionary tracks represents the accretion phase. 
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Fig. 2. Profile of the H abundance as a function of mass for the 
secondary. The abundance profile immediately after accretion is 
shown as solid line. The model with thermohaline mixing is the 
dashed line. The dotted line shows the hydrogen abundance once 
thermohaline mixing has ceased. This model is taken at about 
10 9 years after the initial model. 



The ratio of surface abundances to iron (relative to solar) of 
the model sequences at the end of core hydrogen burning are 
shown in Table Q] Thermohaline mixing has a significant ef- 
fect on the surface abundances. In the model with thermohaline 
mixing all heavy elements are diluted compared to the unmixed 
model, but the mixing does not dilute all the abundances by the 
same amount - the degree of dilution depends on how different 
the abundance of the original material is from that of the ac- 
creted material. For instance, the accreted material is extremely 
carbon-rich and when it is mixed with the pristine stellar mate- 
rial a large change in the abundance results: [C/Fe] is reduced 
from 3.25 in the model without thermohaline mixing to 2.41 in 
the mixed model. Similarly [Na/Fe] and [Mg/Fe] are reduced by 
about 0.8 dex. For other elements the change is less severe. 

3.1. Abundances after first dredge-up 

As a low-mass star begins to ascend the red giant branch, the 
convective envelope deepens and material from the interior of 
the star is drawn to the surface. This is first dredge-up. Material 
that has experienced CN-cycling is brought to the surface and 
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Table 1. Ratio of surface metal abundaces to iron measured relative to solar values. Before accretion, the ratios are zero as we have 
assumed solar-scaled initial values. 



hence the surface carbon abundance drops while the abundance 
of nitrogen increases. If accreted material were to remain un- 
mixed on the surface of the star then during first dredge-up 
it would become mixed with the interior of the star. It is for 
this reason that surveys like that of iLucatello et al.l d2005l) select 
dwarf stars, rather than giant stars. 

Table Q] lists the surface abundances of the two models after 
the occurrence of first dredge-up. For the model without ther- 
mohaline mixing the main effect of the dredge-up is not to bring 
processed material to the surface, but to dilute the formely ac- 
creted material by mixing it with the rest of the envelope - the 
surface abundance of all elements is decreased. The post-dredge- 
up abundances resemble those of the model with thermohaline 
mixing at the end of the main sequence, although the dilution is 
less severe because the dredge-up does not mix as deeply. 

Conversely, in the model that does include thermohaline 
mixing the accreted material had already been mixed during the 
main sequence, so the dredge-up acts like normal, bringing pro- 
cessed material to the surface. Some of the accreted C-rich ma- 
terial was mixed into the interior where it reached temperatures 
high enough for CN-cycling to occur (see the left-hand panel 
of Figure |3). This material is brought to the surface again by 
the dredge-up, causing an increase in [N/Fe] and a decrease in 
[C/Fe]. The [N/Fe] before dredge-up is primarily determined by 
the nitrogen content of the accreted material, whereas [N/Fe] af- 
ter dredge-up is primarily determined by the carbon content of 
the accreted material. 

The thermohaline mixing model shows a slight enhancement 
in [Na/Fe] for a similar reason. The accreted material is rich in 
22 Ne (the abundance of 22 Ne by mass fraction in the accreted 
material is 2.1 x 1CT 3 ). When this material is mixed into the in- 
terior a small amount of this 22 Ne is converted to 23 Na by proton 
capture, even though the reaction rate is small. Since the 23 Na 
abundance is low (its abundance by mass fraction is around 10~ 7 
in the original stellar material) even a small amount of conver- 
sion of 22 Ne contributes significantly to the abundance of 23 Na. 

4. Comparison to observations 

We now compare our models to observations of CEMP stars. 
ILucatello et al.l (120061) published a large sample of CEMP stars 
and their abundances as a function of luminosity. If we assume 
that this data set reflects the chemical evolution of these ob- 
jects, our models should reproduce the trend seen in these data. 
Figure [4] shows the evolution of the [C/Fe] ratio of our mod- 
els against the Lucatello et al. data. The model without ther- 
mohaline mixing (dashed line) displays a large drop in [C/Fe] 
which begins at around log 1() L/L = 0.7 due to the onset of first 
dredge-up, when the accreted material is mixed into the stellar 
interior as the convective envelope deepens. In the model with 
thermohaline mixing (solid line) the accreted material is mixed 
into the interior whilst the star is still on the main sequence, at 
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Fig. 4. Plot of the evolution of the [C/Fe] ratios of the two mod- 
els as a function of luminosity. The solid line is the model with 
thermohaline mixing; the dashed line is the model without ther- 
mohaline mixing. The dotted line is a model with the [C/Fe] 
ratio of the accreted material decreased by one dex and ther- 
mohaline mixing included. Crosses represent th e carbon-rich 
stars of the HERES sample (ILucatello et al. 120061) . The decrease 
in the [C/Fe] ratio of the model with thermohaline mixing at 
log I0 L/L « takes around 10 9 yrs. The excursion to higher lu- 
minosity before [C/Fe] drops corresponds to the accretion phase, 
which is very rapid compared to the rest of the evolution. 



log I0 L/L o « 0. Thus only a small change in [C/Fe] is seen as 
first dredge-up commences. 

Our models are somewhat carbon rich compared to the 
observations. We have chosen a primary that spends a long 
time on the TP-AGB and therefore becomes very carbon-rich. 
Nevertheless, the model with thermohaline mixing is a better 
representation of the data, which show no evidence of a drop 
in [C/Fe] above log 10 L/L = 0.7. The fact that there are no ob- 
served objects with [C/Fe] greater than around 2.5 at luminosi- 
ties less than log I0 L/L Q = 0.7 suggests that accreted material is 
mixed into the companion whilst it is still on the main sequence. 

There is considerable scatter in the Lucatello et al. data. The 
objects in this data set cover a spread of around two dex in metal- 
licity. In the framework of our model, we would expect there 
to be scatter associated with the distribution of the parameters 
of the binary systems (initial primary masses, mass ratios and 
separations). This means there should be a spread in the abun- 
dances of the accreted material. The dotted line in Fig.|4]shows 
the results of a model with thermohaline mixing in which we 
decreased the [C/Fe] ratio of the accreted material by 1 dex. It 
shows the same trend as the original model but at lower [C/Fe], 
covering stars at the lower end of the observed [C/Fe] range. 

Figure |5] shows the evolution of the [N/Fe] ratio for our orig- 
inal model with thermohaline mixing (solid line) against the 
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Fig. 3. Profiles of the 'H, I2 C, 14 N and 16 abundances as a fuction of mass at the end of the main-sequence. The left-hand panel is 
the model with thermohaline mixing and the right-hand panel is the model without thermohaline mixing. The step in the hydrogen 
profile of the thermohaline mixing model is due to the existence of a gradient in the abundance of CNO. The burning rate increases 
rapidly across this region due to the increased abundance in these metals. The models are not quite co-incident in time due to the 
frequency at which the code outputs detailed models. 
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Fig. 5. Plot of the evolution of the [N/Fe] as a function of lumi- 
nosity when thermohaline mixing is included. The solid line is 
for the model computed with the abundances predicted by the 
synthetic code. The dashed line is the model with the [N/Fe] ra- 
tio of the accreted material has been increased by 2 dex. The 
dotted line is the model with the [C/Fe] ratio of the accreted ma- 
terial decreased by 1 dex. T he crosses represent t he carbon-rich 
stars of the HERES sample (iLucatello et al.ll2006l) . 



Lucatello et al. data. When first dredge-up occurs, we see a large 
increase in the [N/Fe] ratio for the reasons outlined above. The 
size of the increase at first dredge-up is dominated by the abun- 
dance of carbon in the accreted material. This model covers the 
upper end of the observed [N/Fe] values but only after dredge- 
up. The observations again show a large spread in [N/Fe] but 
no evidence of a systematic increase at log lo L/L « 0.7. The 
large number of stars with enhanced [N/Fe] at lower luminosi- 
ties in fact strongly suggests that the accreted material was al- 
ready nitrogen-rich. Figure [5] also shows two additional tracks. 
In the model indicated with the dashed line we have boosted 
the [N/Fe] ratio of the accreted material by 2 dex. This track 
passes through the more nitrogen-rich points at low luminosity 
and reaches a similar final [N/Fe] ratio to the original model. 
The second model (dotted line) is the same as shown in Figure|4] 
in which the [C/Fe] ratio of the accreted material is decreased 
by 1 dex. This model reaches a much lower final [N/Fe] ratio 
as should be expected. Hence by varying the carbon and nitro- 



gen abundances of the accreted material the observed range of 
[N/Fe] values in the Lucatello et al. sample can be covered rea- 
sonably well. If thermohaline mixing is not included, the [N/Fe] 
ratio shows a drop at log lo L/L « 0.7, similar to what is seen 
for [C/Fe] in Fig.|4j as the accreted material becomes mixed with 
the convective envelope. 

5. Conclusion 

We have investigated the effects of thermohaline mixing on a 
low-metallicity, low-mass star that has accreted material from 
a carbon-rich companion. This material does not remain on the 
surface of the star as is commonly assumed. A thermohaline in- 
stability leads to the rapid mixing of this material with the rest of 
the star. We find that about 90% of the star experiences mixing, 
while only 0. 1 M at the centre of the star remains unmixed. 

This has a very dramatic effect on the surface abundances of 
the star. The surface abundance of carbon is decreased by nearly 
an order of magnitude relative to the unmixed case. The models 
also display very different abundance patterns after the occur- 
rence of first dredge-up. In the case of the model without thermo- 
haline mixing, the post-first-dredge-up abundances are similar 
to those of the thermohaline mixing model when it has reached 
equilibrium after accretion. On the other hand, the model with 
thermohaline mixing displays a very large increase in its nitro- 
gen abundance after first dredge-up. 

Thermohaline mixing is an important component of models 
of CEMP stars. Without its inclusion, we would expect to see 
an abrupt change in the [C/Fe] ratio as a function of luminosity 
because of first dredge-up mixing the accreted material into the 
stellar interior. Such a change is not observed and thermohaline 
mixing is a possible explanation for this. The observed [N/Fe] 
ratios can be reproduced by either model provided that there is a 
spread in the C and N abundances of the accreted material. 
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